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ABSTRACT 

Numerical simulations of HH jets never show side-entrainment of environmental material into the jet beam. 
This is because the bow shock associated with the jet head pushes the surrounding environment into a dense 
shell, which is never in direct contact with the sides of the jet beam. We present 3D simulations in which a 
side- streaming motion (representing the motion of the outflow source through the surrounding medium) pushes 
the post-bow shock shell into direct contact with the jet beam. This is a possible mechanism for modelling 
well collimated "molecular jets" as an atomic/ionic flow which entrains molecules initially present only in the 
surrounding environment. 

Subject headings: circumstellar matter - hydrodynamics - stars: formation - ISM: HH objects — ISM: jets 
and outflows 



1. INTRODUCTION 

Models of entrainment of molecular, environmental mate- 
rial, in the wings of bow shocks (associated with working sur- 
faces in jets from young stars), are successful in explaining the 
limb-brightened, cavity-like molecular outflows (at least in a 
qualitative way). Analytic (Masson & Chern rn||1993| |Raga 
& Cabrit|[T993] ) and numericaT jRaga et al.||1995[ |Lim et~aL 



2001) models of this so-called "prompt entrainment" scenario 
produce limb-brightened molecular structures that resemble 
the cavity-lile morph ologies observed in objects such as the 
LI 157 outflow (e. g., |Beltran et al.|2004| ). 

However, some outflows from young stars also show high- 
velocity, collimated, jet-like molecular structures. An exam- 
ple of this kind of structure is observed in HH212 (e. g., 
|Codella et al.|2007] ). These jet-like molecular outflows have 
been successfully model ed by assuming that the jet itself is 
initially molecular (e g., |Volker et al.|1999[|Lim et al.|200T 
jMoraghan et al.|2006| ). 



Could these "molecular jets" be the result of environmental 
molecular gas being entrained into an atomic/ionic jet? The 
possibility of having "side entrainment" of molecular mate- 



rial has been studied analytically (Canto & R aga 1991 J and 



numerically (Taylor & Raga| |1995 



Lim et al.||1999l TThese 



models show that if one has a fast, atomic jet beam in direct 
contact (through the sides of the beam) with a molecular envi- 
ronment, a substantial amount of molecular material is indeed 
entrained into the fast flow. 

However, if one computes full simulations of a jet flow, the 
leading bow shock pushes aside the molecular environment, 
so that the sides of the jet beam are never in direct contact 
with the molecular gas. Therefore, the situation nec essary for 
producing side-entrainment of molecular gas (see |Taylor & 
|Raga|1995[|Lim et al.|1999| ) is not obtained. 

In the present paper, we study the possibility of overcom- 
ing this problem by having a low velocity side-motion of the 
environment relative to the jet source. This side-wind could 
represent the motion of the jet source within the surround- 
ing environment. The qualitative effect of the sidewind is de- 
scribed in §2. We have then computed a set of 3D simulations 
of a radiative jet in a sidewind (§3), producing a variety of 

Electronic address: diego.lopez@nucleares.unam.mx 

1 Instituto de Ciencias Nucleares, Universidad Nacional Autonoma de 
Mexico, Ap. 70-543, 04510 D.F., Mexico 



flow morphologies (§4). From the resulting flows we compute 
the amount of environmental material which is pushed by the 
jet flow (§5) and analyze how much material is actually en- 
trained into the jet beam itself (§6). We finally illustrate the 
dependence of our results on the resolution of the numerical 
simulations (§7). The results are summarized in §8. 

2. JET IN A SIDEWIND 

The main problem when trying to incorporate molecular, 
environmental material, into a collimated jet is that the lead- 
ing head of the jet, and possibly also any trailing "internal 
working surfaces", push away the environmental gas into a 
dense shell, which follows the shape of the bow shock wings. 
Because of this, the molecular environmental material never 
reaches contact with the jet beam, and lateral entrainment of 
this material into the jet does not occur. This situation is 
shown in panel a of Figure [T] 





Fig. 1 . — Schematic diagram showing a jet travelling in a stationary envi- 
ronment (top) and in a sidewind travelling parallel to the ordinate (bottom). 
The sidewind pushes the dense, post-bow shock shell into direct contact with 
the jet beam. 

It has been suggested that if the dense shell material were 
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warm enough, it might re-expand into the cavity left by the 
passage of the jet head, and rea ch contact with the jet beam 
(see, e. g., Ra ga & Cabrit|19 93). However, this is f ound not 
to be an important effect in jet simulations (see, e.g., |Raga~et| 
|aLl|1995[ |Lim et aL]|2001| ). Possibly, a stratification of the 
surrounding environment and/or a precession and variabil- 
ity of the jet ejection could l ead to the occurren ce of side- 
entrainment into the jet beam (Cabri fet al.|1997| ). However, 
ntil now the correct combination of parameters for this to oc- 
cur has not been foun d. Other possibil ities have been sug- 
gested. For example, |Lim et aL] ( |2001| ) studied the survival 
of molecules (originally present in the environment) in the 
head of an accelerating jet flow. It is not clear how this high- 
velocity molecular gas in the jet head could end up being en- 
tra ined into the jet bea m. Another possibility was suggested 
by |Raga et aL] ( 2003| ), who proposed that the existence of 
small, dense, moving clumps within the environment might 
be a way of introducing molecular material into the jet beam. 
However, it is not clear that this mechanism leads to molecu- 
lar emission structures that resemble the observations. 

The possibility that we study in this paper is that the pres- 
ence of a side- streaming environment pushes the bow shock 
wing (and the post-bow shock, dense shell) against the jet 
beam, as shown in panel b. of Figure [T] The side- streaming 
could be the result of the motion of the jet source through 
the molecular cloud, and w ould have velocities of at most a 
few km s -1 . Both anal ytic ([Canto & Raga||l995[|Raga et al. 
2009a ) and numerical (|Lim & Ra ga 1998; Masciadri & Raga 



TABLE l 
Model characteristics. 
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Ciardi 



et al.||2008| ) models of jets in sidewinds have 
been computed previously. These models have mostly been 
applie d to jets within expanding H II regions ( [Ciardi et al.| 
2008) or to jets embedded in an isotropic stellar wind (Raga 
et al.|2009a| . The relatively high sidewind velocities relevant 
for these cases (~ 10-30 km s -1 for an expanding H II region 
and up to ~ 1000 km s -1 for a stellar wind), can produce jets 
with strongly curved jet beams. 

For lower sidewind velocities (not studied in the papers 
cited above), only a weak curvature will be produced as a re- 
sult of the jet/sidewind interaction. However, the bow shock 
wing will still be pushed against the body of the jet, as shown 
in the schematic diagram of Figure [T] We focus on this 
regime, in which the dense shell of swept up environmental 
material is pushed into contact with the body of the jet beam 
(therefore allowing the entrainment of environmental gas into 
the beam), but a relatively straight jet path is still obtained. 
The numerical simulations which we have carried out are de- 
scribed in the next section. 

3. THE NUMERICAL SIMULATIONS 

We have computed a set of 3D gasdynamic simulations of 
jet/sidewind interactions. All of them have been computed in 
a (2,0.5,0.5) x 10 17 cm cartesian grid. The jet is injected at 
x = (in the centre of the boundary plane of the computational 
grid), with a velocity parallel to the x-axis. A sidewind is 
injected in the y = plane, with a velocity directed along the 
v-axis. A reflection boundary condition is applied on the x = 
boundary outside the jet beam, and transmission conditions 
are applied on all of the other boundaries except the y = 
plane (in which the sidewind is injected). 

An initially neutral, top-hat jet of velocity V/, density rij, 
radius r 7 = 10 15 cm and temperature 7} = 10 3 K moves into 
an initially uniform, neutral environment with a density n a = 
200 cm -3 , temperature T a = 10 K and sidestreaming velocity 
v a . A set of models with different values of Vj, rij and v a has 
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been computed, with the parameters given in Table 1 . 

The simulations were carried out with the "Yguazu-a" 
code (Ra ga et aT||2000| ), solving the 3D gasdynamic equa- 
tions together with a continuity /rate equation fo r neutral H. 
The par ametrized cooling function described by Raga et al. 
(2009b) is included in the energy equation. Also, we inte- 
grate an equation for a normalized passive scalar g with which 
we distinguish between the ambient and jet medium. If the 
scalar was positive it indicated that the material was initially 
medium material, while if it was negative it was jet material. 
For example, if we had only ambient medium material g = 1 , 
or if we had pure jet material g = -1. With the use of this scalar 
we were able to calculate the amount of mixing between the 
ambient medium and the jet material. For this, we defined the 
mixing mass fraction as f m = (g+ 1)/2. It is clear that for pure 
ambient medium f m = 1 ; while for pure jet material f m = 0, in- 
termediate values of f m indicates that there was material from 
both the ambient medium and the jet mixed together. The 
case in which we have 99% ambient medium mixed with 1% 
of jet material (f m = 0.99), is defined as "99% mixing frac- 
tion"; while the 1% ambient medium case (f m = 0.01) as "1% 
mixing fraction". 

A 6-level binary adaptive grid has been used with three 
different maximum resolutions Ax = 3.91, 1.95 and 0.98 x 
10 13 cm (along the three axes). We have called these the 
"low", "medium" and "high" resolutions (labeled in Table 1 
with letters lr, mr and hr, respectively). 



4. THE RESULTING FLOW STRATIFICATIONS 

As an example of the flows resulting from our simulations, 
in Figure[2]we show xy-cuts showing the time-evolution of the 
mid-plane density stratification obtained from model al (see 
Table 1). This figure also shows 2 contours, corresponding to 
the 99% mixing fraction (f m = 0.99, outer contour), and 1% 
mixing fraction (f m = 0.01, inner contour). 

In model al, the xy-cuts show a side-to-side asymmetry that 
is a direct result of the fact that the environment is flowing 
along the +j-direction. This asymmetry is seen as a distortion 
of the leading bow shock, and as a penetration of environmen- 
tal material to regions close to the jet beam in the up-sidewind 
(-y) direction. 

In Figure [3] we show single time frames of the xy-midplane 
density stratifications obtained from the 9 models of Table 1 . 
In each of the three columns of Figure [3] (from top to bot- 
tom), we see the effect of increasing velocity of the sidewind 
(the top, centre and bottom frames correspond to v a = 2, 5 
and 10 km s -1 , respectively, see Table 1). It is clear that 
for higher values of v a , the environmental material penetrates 
more strongly towards the jet beam in the up-sidewind region. 
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Fig. 2. — Density stratifications (color scale, given in g cm -3 by the top 
bar) and mixing fractions (contours) obtained from model al (see Table 1). 
The side-streaming environment flows along the vertical axis of the plots. 
The displayed stratifications correspond to cuts on a plane that includes the 
outflow axis and the sidewind, at integration times t = 100 (top), 200 and 
300 yr (bottom). The two contours correspond to the 99 % mixing fraction 
(fm = 0.99, outer contour); and the 1 % mixing fraction case (f m = 0.01, inner 
contour). The axes are labeled in units of 10 16 cm. 

A comparison of the first and second columns of Fig- 
ure [3] shows that if one increases the jet velocity from vj = 
150 km s -1 (models al, bl and cl) to 300 km s -1 (models 
a2, b2 and c2, see Table 1) the resulting density stratifications 
and mixing fractions remain qualitatively unchanged, show- 
ing slightly more pronounced asymmetries for the higher jet 
velocity. In particular, it is clear that in regions close to the 
source the post-bow shock dense shell touches the jet beam in 
all of the v a = 10 km s -1 models (bottom row of Figure (3J). 

Finally, we also see that if one increases the jet density 
(from rij = 1000 cm -3 for the models in the first column to 
5000 cm -3 for the models in the third column, see Table 
1), less penetration of the environmental material into the jet 
beam is obtained. This effect can be seen as broader 99% 
pure jet material regions (limited by the inner contour) in the 
a3, b3, c3 models (compared to al, bl and cl). 

5. ENTRAINED MATERIAL 

As we have described in ^3] from our simulations we ob- 
tain the environmental to total mass mixing fraction f m as a 
function of position and time. Using this mixing fraction, we 
compute the jet mass loss rate associated with motions along 
the x-axis : 

Mj(x) = J J [l-f m ]pudydz, (1) 
and the mass rate of the entrained material 

M AM (x) = J J f m pudydz, (2) 



where f m , p and u are the 3D mixing fraction, density and x- 
velocity (respectively), obtained for a given integration time 
t. 

The mass loss rates obtained in this way for models al and 
cl are shown in Figure [4] and [5] (for times t = 50 to 300 yr). It 
must be noted that calculating the position of the jet head was 
not trivial, it required fine-tuning, was not obtained exactly 
the same for each model, and the environmental mass rate 
from the entrained material (Mam) was extremely sensitive to 
it. Thus, in order to be consistent in our analysis, we exclude 
the head of the jet from our discussion. For model al (which 
has a v a = 2 km s -1 sidewind, see Table 1), we see that, at 
all times Mam monotonically grows along the jet axis, and 
that a maximum is reached at the position of the jets head 
(top panel of Figure [4]). For this model, the Mam/Mj ratio 
also grows with distance x from the source, having values of 
~ 5 — )> 8 x 10~ 3 at the middle of the length of the jet (at a given 
integration time), and reaching values of ~ 1 at the head of the 
jet (bottom panel of Figure]?]). 

For model cl (which has the same parameters as model al 
except for a v a = 10 km s -1 sidewind, see Table 1), a qualita- 
tively similar behavior is obtained for the ambient mass rate, 
but with higher values of Mam at all times and positions along 
the jet (top panel of Figure pj). The M A m/Mj ratio has values 
of ~ 1.1 1.3 x 10" 2 at the middle of the length of the jet 
(for all times), a factor of ~ 2 larger than the values obtained 
for model al (also except for the head of the jet, which we 
excluded from the analysis due to lack of consistency). 

We now calculate the mean velocities associated with the 
jet and environmental mass rates : 



vj(x) = 



1 



Mj(x) 



vam(x) = 



J J [l-/m] pu 2 dydz, 
J J f m pu 2 dydz. 



(3) 



(4) 



M AM (x) 

where Mj(x) and Mam are given by equations and ([2]). 

The mean velocities obtained in this way for models al and 
cl are shown in Figure [6] We see that the mean forward ve- 
locity Vj of the jet material is of 150 km s -1 at x = (correctly 
coinciding with the injection velocity, see Table 1). For model 
al (top panel of Figure [6]), the jet velocity drops close to the 
position of the head to ~ 120 km s -1 . The drop in Vj occurs 
closer to the jet source in model cl (bottom panel of Figure|6]). 
In both models (al and cl), the mean forward velocity Vam 
of the environmental material grows along the length of the 
outflow, with a velocity of « 120 km s -1 at the jet head. At 
all times, we also see a peak in Vam at x < 10 16 cm (the po- 
sition of this peak being at larger distances from the source 
for longer integration times). This peak is associated with a 
region (close to the jet source) in which the post-bow shock 
shell touches the jet beam since early evolutionary times (see 
Figure [2]). 

In Figure [7] we show the environmental mass rate as a 
function of distance from the source obtained from all of 
the models of Table 1 (only the results for a single integra- 
tion time are shown for each model). All of the models give 
Mam ~ 10 17 g s -1 close to the head of the jet. Lower mass loss 
rates are obtained closer to the outflow source. 

The models with Vj = 300 km s -1 (a2, b2 and c2, cen- 
tral panel) show basically identical M A m(x) dependencies for 
x < 5 x 10 16 cm, and a clear spike at x sp w 10 16 cm, but as 
we shall see corresponds to environmental material which is 
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Fig. 3. — Density stratifications (color scale, given in g cm -3 by the top bar) and mixing fractions (contours) obtained for the models listed in Table 1. The 
displayed stratifications correspond to cuts on a plane that includes the outflow axis and the sidewind, at integration times t = 200 yr (for models al, bl, cl, a3, 
b3 and c3) or t = 100 yr (for models a2, b2, c2). The two contours correspond to the 99 % mixing fraction (f m = 0.99, outer contour); and the 1 % mixing fraction 
case (f m = 0.01, inner contour). The axes are labeled in units of 10 16 cm. 
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Fig. 4. — Top panel: Mass rates associated with the jet (Mj, black lines) 
and environmental material (Mam, coloured lines) as a function of position 
x along the outflow axis, obtained from the flow stratifications of model al 
(see Table 1) at integration times t\ = 50 (red lines), ti - 100 (green), t$ = 150 
(blue), t 4 = 200 (pink), t 5 = 250 (cyan) and t 6 = 300 yr (yellow). Bottom 
panel: the Mam/Mj ratio as a function of x for the same integration times. 



not entrained primarily by the sidewind, but rather from the 
leading bow-shock, thus is not of interest for this study. The 
Vj = 150 km s -1 models (top and bottom panels) show a re- 
gion close to the outflow source (with x < 2 x 10 16 cm) in 
which larger Mam are obtained for larger sidewind velocities 
{y a = 2 km s -1 for models al and a3, 5 km s -1 for bl and b3 
and 10 km s -1 for cl and c3, see Table 1). As we will see in the 
following section, in this region we are seeing environmental 
material that is directly entrained into the jet beam (due to the 
presence of the sidewind). 

Figure M shows the environmental mean velocity (vam, see 
equation as a function of position for all of the computed 
models (see Table 1). In the Vj = 150 km s -1 models (upper 
and lower panels of Figure [8]), regardless of the sidewind ve- 
locity the correspondent vam grows from low (~ 1 km s -1 ) 
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Fig. 5. — The same as Figure 4, but for model cl (see Table 1), 
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close to the source, up to values comparable to the jet velocity 
(~ 250 km s -1 ) at the head of the jet. 

Very similar results are obtained for the models with Vj = 
300 km s -1 (central panel), except for the presence of a spike 
at x sp « 10 16 cm. Since the corresponding vam &tx sp for these 
models (~ 0.03v 7 ) was lower than that from the models with 
vj = 150 km s -1 (where vam ~ 0.04-0.20v 7 ), we had an in- 
sight that the material at the spike of Figure [7] (for models 
with vj = 300 km s -1 ) corresponded to slowly moving mate- 
rial which was mostly entrained by the leading bow- shock and 
not due to the presence of the sidewind. 

Also in Figure [Sj there is evidence of "high velocity materi- 
alO. Our models with Vj = 150 km s -1 showed the presence of 
a fast region close to the source (x < 2 x 10 16 cm). In this re- 
gion vam has values of ~ 40 — » 70 km s -1 (of order ~ 50% of 
the Vj). These large vam values close to the source are associ- 
ated with the direct entrainment of the dense post-bow shock 
shell into the jet beam due to the sidewind. 
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Fig. 6. — Top panel: Average velocity associated with the jet (vy, black 
lines) and the environmental material (vam, coloured lines) as a function of 
position x along the outflow axis, obtained from the flow stratifications of 
model al (see Table 1) at integration times t\ = 50 (red lines), £2 = 100 (green), 
t 3 = 150 (blue), t 4 = 200 (pink), t 5 = 250 (cyan) and t 6 = 300 yr (yellow). 
Bottom panel: The same as in the top panel, but for model cl (see Table 1). 
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Fig. 7. — Environmental mass rate Mam (see equation^ computed for all 
models (see Table 1) for an integration time of 200 yr for models al, bl, cl, 
a3, b3 and c3 and of 100 yr for models a2, b2 and c2. 
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Fig. 8. — Environmental average velocity (vam) (see equation|4j computed 
for all models (see Table 1) for an integration time of 200 yr for models al, 
bl, cl, a3, b3 and c3 and of 100 yr for models a2, b2 and c2. 



6. FAST ENTRAINED MATERIAL 



In order to study the environmental material which is di- 
rectly entrained into the jet beam, we compute the mass rate 
and average velocity of the environmental material that moves 
along the jet axis at velocities >Vj/2 (in other words, with ve- 
locities larger than 75 km s -1 for models al, bl, cl, a3, b3, c3 
and larger than 150 km s -1 for models a2, b2 and c2). We call 
these the "high velocity" mass rate M A M,h and average veloc- 
ity V A M,h- 

The values of M AM ^(x) obtained for all models are shown 
in Figure [9] We see that in the region close to the outflow 
source (x < 2 x 10 16 cm) the low velocity jet models (with 
vj = 150 km s" 1 , models al, bl, cl, a3, b3, c3) have environ- 
mental mass rates which monotonically grow with increas- 
ing values of the sidewind velocity (y a ). For v a = 10 km s -1 
(models cl and c3), the mass rate in this region has values of 
M A M,h ~ 5 x 10 14 g s" 1 , of the order of ~ 0.5% of the mass 
loss rate of the jet. Lower values of M A M,h are obtained for 
the v a = 5 (bl and b3) and v a = 2 km s -1 models (al and a3). 
Finally, we see that the values of M A M,h in the region close to 
the source are much lower (M A M,h ~5x 10 13 g s -1 ) for the 
vj = 300 km s -1 models (a2, b2, c2, central panel of Figurej9J). 
This confirms the fact that the material from the spike in Fig- 
ure [7] (for models with vj = 150 km s -1 ), corresponds to envi- 
ronmental material which was not entrained in its majority by 
the sidewind (and so, is not of interest for this study). 
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Fig. 9. — Mass rate of the environmental material which moves with "high" 
velocities along the jet axis (M A M,h)- For model al, bl, cl, a3, b3, c3 we show 
the material that is moving with v A m >75 km s _1 , and with vam > 150 km s _1 
for models a2, b2 and c2. 

All of the models show a strongly increasing M AM ^(x) for 
x > 2 x 10 16 cm. This growth is associated with the fact that 
the motion of the post-bow shock shell becomes progressively 
more forward directed as we approach the head of the jet. 

Therefore, there are two components of the fast (with axial 
velocities > Vj/2) ambient material : 

• material originating in the region in which the post-bow 
shock shell is in contact with the jet beam, 

• material in the post-bow shock shell in the region close 
to the jet head. 

The spatial distribution of these two components is shown in 
Figure [10| with the mid-plane spatial distribution of the den- 
sity pf m of the fast (y x > Vj/2) ambient material obtained from 
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models al and cl (for an integration time of 300 yr). These 
stratifications show that the fast ambient material is confined 
to a region within or in near contact to the jet beam. 
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Fig. 10. — Environmental mass fraction times the density stratification 
(p X f m , colour scale, given in g cm -3 by the top bar), of the fast ambient 
material obtained from models al (top panel) and cl (bottom panel), for an 
integration time of 300 yr. The axes are labeled in units of 10 16 cm. 



7. A RESOLUTION STUDY 

In order to illustrate the effect of the numerical resolution, 
we have computed two of the models (al and cl, see Table 
1) at three resolutions: Ax = 3.91, 1.95 and 0.98 x 10 13 cm 
(along the three axes). In Figure [TTJ we show the density 
stratifications obtained from model al at these three resolu- 
tions, for a t = 200 yr integration time. 

It is clear that more complex structures are obtained for in- 
creasing resolutions. However, the result that the post-bow 
shock shell is swept into contact with the jet beam (as a result 
of the side- streaming environment) is present at all resolutions 
(see Figure [TT]). Therefore, side-entrainment into the jet beam 
occurs regardless of the resolution of the simulations. 

we show the total environmental mass rate 



In Figure 12 



(Mam, top frame) and the mass rate of the fast entrained mate- 
rial (M AM:h , bottom frame) as a function of distance from the 
source, computed from the t = 200 yr stratification of model 
al . The values of Mam and Mam,h are similar at all resolutions 
in the region close to the jet head. 

On the other hand, in the region close to the source (where 
the dense, post-bow shock shell is being entrained into the jet 
beam), the environmental mass rate shows a stronger depen- 
dence on the resolution. In particular, we have a factor of ~ 10 
decrease in Mam,h for a factor of 4 increase in the resolution of 
the simulation. This result is to be expected (given the lower 
numerical diffusion of the higher resolution simulations), and 
its implications are discussed in the following section. 

It is clear that in this resolution study a relatively small 
range (of a single octave) of resolutions is explored. This is 
a result of the fact that simulations with less than ~ 10 grid 



Fig. 11. — Density stratifications (color scale, given in g cm -3 by the 
top bar) and mixing fractions (contours) obtained from model al at three 
different resolutions for an integration time of 200 yr (see Table 1). The two 
contours correspond to the 99 % mixing fraction (f m = 0.99, outer contour); 
and the 1 % mixing fraction case (f m = 0.01, inner contour). The axes are 
labeled in units of 10 16 cm. The three panels correspond to "low" (case in 
which Ax = 3.91 x 10 13 cm along the three axes, top panel), "medium" (Ax = 
1.95 x 10 13 cm), and "high" (Ax = 0.98 x 10 13 cm, bottom panel) resolutions. 
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Fig. 12. — Total environmental mass rate (Mam, top panel) and the mass 
rate of the fast entrained material (M AM ^, bottom panel), as a function of 
distance from the source, computed for the three different resolutions from 
model al (at an integration time t = 200 yr). 



points across the jet beam are basically meaningless, and that 
our highest resolution simulation (model hr) resolves the jet 
beam diameter with ~ 100 points (at the highest resolution of 
the adaptive grid). This resolution is "competitive" in terms 
of recent 3D astrophysical jet simulations, for example Rossi 
et al. ( 2008 ) computed 3D simulations of entrainment in rel- 
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ativistic jets resolving the jet diameter with ~ 20 points). In 
the near future we will be able to carry out 3D simulations 
with higher resolutions, hopefully reaching ~ 300-500 points 
across the jet diameter. This will extend the range in resolu- 
tions to ~ 1.5 octaves (which is still not very impressive). 

The situation is therefore somewhat hopeless. If we take 
the number of grid points across the jet diameter as an esti- 
mate of the Reynolds number of the simulation, we see that at 
the present time we can only expect to reach Re~ 500. Such 
Reynolds numbers are ~ 2 orders of magnitude below the 
Reynolds number (of ~ 50000) necessary to reach the "high 
Re" regime, in which the entrainment rate becomes indepen- 
dent of Re ( [Birch & Eggers 1972). Therefore, it is unclear 
if this "high Re" regime (relevant for astrophysical jets) will 
ever be reached by 3D gasdynamic simulations. Regardless of 
this fact, exploring the development of the Kelvin-Helmholtz 
instabilities (that initiate the turbulent entrainment) at increas- 
ing spatial resolutions is probably still a worthwhile pursuit. 
Different efforts in this direction have been done by Micono 
^nr] ( |2000a|bl ); |Xu et al.| ( |2000| ); |Rosen & Hardee| ( |2002| ). 

An eventually more fruitful approach might be to incorpo- 
rate a "turbulence parametrization scheme", which has to be 
calibrated with laboratory experiments in order to provide the 
correct mass entrainment rate. The application to astrophys- 
ical jets of a simp le, "a-model" turb ulence parametrization 
was discussed by Canto & Raga (1991). A more evolved, "k- 
e" scheme was explored by [Falle (1994]). 

The problem of this approach of course is that the turbu- 
lence parametrization schemes are normally calibrated with 
non-radiative, laboratory jet flow experiments. Therefore, it 
is not clear whether or not the resulting parametrization is ap- 
propriate for modeling the entrainment in radiative jets. A 
possible solution to this problem would be to attempt to cal- 
ibrate turbulence parametrization schemes with newly avail- 
able experiments of radiative plasma jets (see, e. g. Ample- 
ford et aL||2008| ). The relevant experimental data, however, 
are still not available. 

8. CONCLUSIONS 

It has been suggested that the molecular emission observed 
along well collimated, jet-like outflows from young stars 
might be the result of entrainme nt of molecula r environmen- 
tal material into the jet beam (Canto & Raga 1991). How- 
ever, this side-entrainment has never been seen in numerical 
simulations of HH jets, because the environmental material 
is pushed into a dense, post-bow shock shell which does not 
touch the jet beam. 

In this paper, we present 3D numerical simulations of a 



jet in a sidewind, with sidewind velocities in the v a = 2 — >• 
10 km s -1 range (the lower part of this range being consistent 
with the peculiar motions of T Tauri stars). We find that for 
many parameter combinations the sidewind pushes the post- 
bow shock shell into direct contact with the jet beam (see Fig- 
ure [2] and [3]). In this region of contact, side-entrainment of 
environmental material into the jet beam does take place. 

In our simulations, the side-entrainment results in mass 
rates Mam.h ~ 5 x 10 14 g s -1 , corresponding to ~ 0.5 % of 
the mass loss rate Mj of the jet (see 95b. If the molecular, 
environmental gas is not dissociated during the process of 
side-entrainment, this would result in a molecular fraction of 
~ 0.5 % within the jet beam, which would result in molecular 
column densit ies high enough to produce observable molecu- 
lar emission (Raga et al. 2005]). 

Our present simulations do not include the chemistry of the 
entrained material, so that we are not able to see whether or 
not the molecules in the side-entrained material actually sur- 
vive the entrainment process. However, the fact that the side- 
entrained material has been shocked by the slow-moving far 
bow shock wings, and that the region of contact between the 
shocked environment and the jet remains cool (at tempera- 
tures of ~ 10 3 K) indicates that molecules indeed might be 
entrained into the jet beam without being dissociated. 

Furthermore, our simulations do not describe correctly the 
entrainment of the post-bow shock shell into the jet beam. 
An indication of this is the fact that we obtain entrained mass 
rates that strongly depend on the spatial resolution of our sim- 
ulations (see In order to overcome this problem, it will 
be necessary to go to much higher resolutions, in order to re- 
solve the Kelvin-Hel mholtz instabilities that produce the side- 
entrainment (see e. g |Micono et al.|2000a|b|), or to use a "tur- 
bulence parametrization recipe" dCanto & Ragall 19911 Falle 

Even though our simulations do not fully describe the 
side-entrainment process, they conclusively show that a 
side- streaming environment (reflecting the motion of the jet 
source) will push the post-bow shock shell into direct con- 
tact with the jet beam. This then provides the conditions in 
which molecular material will be entrained into the fast jet 
beam, giving "parmetrized turbulenc e" jet models -initially 
suggested for astrophysical jets by [Bicknellj ( 1984| -, a new 
life as a possible explanation for the molecular jets observed 
in star forming regions. 
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